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The analysis of deeply virtual meson production is extended to neutrino-production of the pseudo- 
Goldstone mesons (ji, K, rj) on nucleons, with the flavor content of the recoil baryon either remaining 
intact, or changing to a hyperon from the SU(3) octet. We rely on the SU(3) relations and express 
CN| ' all the cross-sections in terms of the proton generalized parton distributions (GPDs). The corre- 

sponding amplitudes are calculated at the leading twist level and in the leading order in a s , using 
a phenomenological model of GPDs. We provide a computational code, which can be used for eval- 
(N- uation of the cross-sections employing various GPD models. We conclude that these processes can 

, be studied in the experiment Minerva at FERMILAB, which could supplement the measurements 

(£) ■ at JLAB helping to extract the GPD flavor structure from data. 
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INTRODUCTION 



During the last decade the notion of generalized parton distributions (GPDs) became a standard theoretical tool to 
'qj describe the nonperturbative structure of the hadronic target. These new objects, being an special case of the general 
i <~| Wigner distributions, contain rich information about the nonperturbative dynamics of the target structure, such as form 
factors, ordinary parton distribution functions (PDF), fractions of the spin carried by each parton, etc. (see e.g. recent 
reviews in [H-H). In hard exclusive reactions in the Bjorken kinematics, due to collinear factorization @, II] the amplitude 
j>! ■ of many processes may be represented as a convolution of the process-dependent perturbative coefficient functions with 
| target-dependent GPDs. While a mo del- independent deconvolution and extraction of GPDs from data is in general 
1 impossible 1 , nevertheless data help to constrain the available models of GPDs. 

Currently the main source of information about GPDs are the electron-proton measurements done at JLAB and HERA, 
'^T . in particular deeply virtual Compton scattering (DVCS) and deeply virtual meson production (DVMP) |3l la. M. la. I12h23| . 
C^> ' A planned CLAS12 upgrade at JLAB will help to improve our understanding of the GPDs |23j |. 

Having only data on DVCS one cannot single out the flavor structure of GPDs. The process of DVMP potentially 
^1 ■ is able to disentangle the flavor structure of the GPDs, since different mesons are sensitive to different GPD flavor 
combinations 0, However, the practical realization of this program suffers from large uncertainties. In the HERA 
kinematics {xsj < 10~ 2 ), one is close to the saturation regime, where gluons dominate, and as was discussed in [26j |. NLO 
corrections in this kinematical range become large due to BFKL-type logarithms. As a consequence, for the description 
of exclusive processes one should use models where the saturation is built-in [27rl30l |. At the same time, in the JLAB 



X 



1 kinematics, the range of Q 2 is quite restricted, and the GPDs extracted from DVCS may be essentially contaminated by 
higher-twist effects. In case of the DVMP, an additional uncertainty comes from the distribution amplitudes (DA) of the 
produced mesons: while there is a lot of models, only the DAs for it and rj were confronted with data [31rl33j (see also the 
recent review in [UHil). For heavier mesons (p, ui,(f)), the DAs are completely unknown, because their partonic structure 
is controlled by confinement, rather than by the chiral symmetry as for Goldstone mesons. In the general case it is not 



1 An exception is the process of double deeply virtual Compton scattering discussed in [9l 4lll , however vanishingly small cross-sections make 
it unreachable at modern accelerators. 
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even known if the corresponding DAs should vanish at the endpoints of the fractional light-cone momentum distribution, 
as is required by collinear factorization, or if the amplitude gets a contribution from transverse degrees of freedom. 

From this point of view, consistency checks of GPD extraction from JLAB data, especially of their flavor structure, 
are important. Neutrino experiments present a powerful tool, which could be used for this purpose. Up to recently, the 
high-precision exclusive neutrino-hadron differential cross-sections were available only in the low-energy region, where the 
physics is described by s-channel resonances [36-40]. In the high-energy regime, due to the smallness of the cross-sections 
and vanishingly small luminosities in the tails of the neutrino spectra, the experimental data have been available so far only 
either for inclusive or for integrated (total) exclusive cross-sections. The situation is going to change next year, when the 
high-intensity NuMI beam at Fermilab will switch on to the so-cal led middle-energy (ME) regime with an average neutrino 
energy of about 6 GeV. In this setup the Minerva experiment |4lT ] should be able to probe the quark flavor structure 
of the targets. Potentially, NuMI neutrino beam may reach energies up to 20 GeV, without essential loss of luminosity. 
Even higher luminosities in multi-GeV regime can be achieved at the planned Muon Collider/Neutrino Factory [i^l - liij ] . 

One can access the GPD flavor structure in neutrino interactions by studying the same processes as in ep collisions and 
employing the difference of the weak and electromagnetic couplings for the vector current. An example of such a process 
is the weak DVCS discussed in 45]. However, the weak DVCS alone is not sufficient to constrain the flavor structure. 
Moreover, the typical magnitude of the cross-sections for such processes is tiny, of the order ~ 10~ 44 cm 2 / GeV 4 . 

The i/DVMP measurements with neutrino and antineutrino beams are complementary to the electromagnetic DVCS. 
In the axial channel, due to the chiral symmetry breaking we have an octet of pseudo-Goldstone bosons which act as 
a natural probe for the flavor content. Due to the V — A structure of the charged current, in v DVMP one can access 
simultaneously the unpolarized GPDs, H, E, and the helicity flip GPDs, H and E. Besides, important information on 
flavor structure can be obtained by studying the transitional GPDs in the processes with nucleon to hyperon transitions. 
As was discussed in (46{, due to SU(3) flavor symmetry, these GPDs can be related to the ordinary diagonal GPDs in the 
proton. 

The paper is organized as follows. In Section [TT] we evaluate the Goldstone meson production by neutrinos on nucleon 
targets. The main result of this section is Table Q] and Eqns I6HU1 In Section IIII1 for the sake of completeness we sketch 
the properties of the GPD parametrization which will be used for evaluations. In Section ITVl we present numerical results 
and make conclusions. 

II. CROSS-SECTION OF THE i/DVMP PROCESS 

The DVMP process in the vector channel has been studied in [24], [25| . In the leading order in a s , the hard coefficient 
function gets contributions from the diagrams shown in Figure [TJ 




FIG. 1: Leading-order and leading twist contributions to the DVMP hard coefficient functions. 
Evaluation of these diagrams is straightforward and yields for the amplitude of the process, 

Tm = (/ dz^^HlNMTN^ (1) 

where N(p), N(p) are the spinors of the initial/final state baryon, 4>m{z) is the normalized to unity distribution amplitude 
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of the produced meson, Jm is the decay constant of the corresponding meson, Hj^N (pa) TN (pi) is a symbolic notation 
for summation of all leading twist GPDs contributions (defined below) , and T-L T M are the convolutions of the GPDs Hy of the 
target with the proper coefficient function. Currently, the amplitude of the DVMP is known up to NLO accuracy [47], |48| . 
Extension of the analysis of [24], [H| to neutrinos is straightforward. In contrast to electroproduction, due to V — A 
structure, the amplitudes acquire contributions from both the unpolarized and helicity flip GPDs. 
In the leading twist, four GPDs, H, E, H and E contribute to this process. They are defined as 



2^ 



dze ixp+z (B> 
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where P = P i + p 2 , A = P2 — p\ and £ = — A+/2P+ ps XBj/{2 — XBj) (see e.g. [13] for details of kinematics). In the 
general case, when A ^ B, in the right-hand side (r.h.s.) of Eqs. J5|, © there might be extra structures which are 
forbidden by T-parity in the case of A = B 20]. In what follows we assume that the initial state A is either a proton 
or a neutron, and B belongs to the same lowest SU(3) octet of baryons. In this case, all such terms are parametrically 
suppressed by the current quark mass m q and vanish in the limit of exact SU(3), so we will disregard them. Since in 
neutrino experiments the target cannot be polarized due to its large size, it makes no sense to discuss contributions of 
the transversity GPDs Ht, Et, Ht, Et- Also, in this paper we ignore the contributions of gluons, because in the current 
and forthcoming neutrino experiments the region of small xsj <C 1 is not achievable, so the amplitude (QJ simplifies to 



8ni a s f M 



dz 



4>m{z) 



UmN (pa) 7+75^ (pi) + J^EmN (p 2 ) 75 iV (pi) 



2tojv 



(4) 



^M^(p 2 )7+^(Pi) 
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In table U the corresponding amplitudes are listed for each final state M. It is restricted to the cases of either protons 
or neutrons in the initial state, and only baryons from the lowest lying octet in the final state. We used ordinary SU (3) 



relations 146] to relate the transitional GPDs (Y 



p) to the proton GPDs (p 







As was mentioned in I4£ 



these relations for the GPD E can be broken due to the different masses of pion and kaon in the i-channel. Also, the 
SU(3) relations can be inaccurate at sm&M-x Bj (high energy), due to different intercepts of the ir/p and K/K* Regge 
trajectories 2 0, Besides, as was discussed in [26|, in the small- a; regime NLO corrections become large due to 
BFKL-type logarithms, and a lot of care is needed to make a systematic resummation and avoid double counting. For 
this reason, in what follows we restrict our consideration to the moderate energy range (xbj ^ 0.1). For a neutron target, 
in the right panel of table [J we flipped H u / n —> H^/p, Hd/ n — > H u / p , so all the GPDs are given for a proton target. The 



corresponding constants <?/ should be understood as neutral current couplings g A and g J v for the DVMP form factors H 



J 



£, and H, £ respectively. Electroproduction data |2J,l25| correspond to the mere change of the charges, gA 



0, 9 f v 



2 Similar results may be obtained in the framework of the dipole model [27l430l |. which is valid for very small x < 10 -2 : The amplitude gets 
a substantial contribution from the endpoint region with a or a ~ m 2 /Q 2 . This asymmetry in the a contribution depends strongly on the 
quark mass and obviously breaks the SU (3) symmetry. 
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TABLE I: List of the DVMP amplitudes Hm, Em, Hm, Em for different final states. For a neutron target, in the r.h.s. we flipped 
H u / n — > H d / p , H d / n — > H u /„, so all the GPDs are given for a proton target. To get £ , H, E, replace H with E, H, E respectively. 
The corresponding constants gj should be understood as neutral current couplings g^ A for the DVMP form factors H, £, and as g(r 
for Iri, E. Vij are the CKM matrix elements. c± is a shorthand notation c±(x, £) — l/(a; ± £ =p iO) for the leading order coefficient 
function. NLO corrections to the coefficient functions may be found in [47], |4g] . For the sake of brevity, we did not show the 
arguments (a;, £, t, Q) for all GPDs and omitted the integral over the quark light-cone fraction J dx everywhere. 
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While in the table we listed, for the purpose of reference, all 41 amplitudes, only 12 of them are independent due 
to the SU(3) symmetry 3 . This agrees with the fact that all the amplitudes are linear combinations of 6 functions 



3 This follows from number of irreps in 8x8=l+8+8+10+10*+27 for a given J p 
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J dx H Ut d, s (x, £, t) c±(x, £) and 6 functions J dx H u ,d,s {x, £, t) c±(x, £) for the axial and vector channels respectively 4 . This 
implies a large number of relations between different cross-sections, some of which are obvious. For example, comparing 
different elements of the Table |H we may get 5 : 
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Other relations can be extracted using the identity 

\A + B\ 2 + \A- B\ 2 = 2(\A\ 
For example, fixing 

A=(H U -H d )c, B — (H u — H d ) c + , 
arrive at the relation between the Cabibbo-suppressed and Cabibbo-allowed cross-sections 



we 
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The corresponding neutrino cross-sections for charged and neutral currents read, 
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In analogy to the electro- and photoproduction processes, it makes sense to introduce the cross-section of the subpr 
W/Z + p —> M + p, which has the form, 
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4 As one can see from [lOl the GPDs E, E always contribute in bilinear combinations H*H — E* E with the same coefficients as the GPDs 
H, H , so this does not change the total count of independent cross-sections 

5 For example, the cross-sections of up — > /i~tt + p and i/n -> fj, + w~n are equal because the amplitudes of subprocesses W + p — > ir+p and 
W~n — > ir~n coincide due to isospin symmetry (/ = 3/2 state). 
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Iii neutrino experiments the target is unpolarized, so |T A /| can be simplified to 
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III. GPD AND DA PARAMETRIZATIONS 



As was mentioned in the introduction, in the case of the DVMP a large part of the uncertainty comes from the DAs of 
the produced mesons. In spite of many mo del- dependent estimates, so far DAs are poorly known. Experimentally, only the 
distribution amplitudes of the pions and ?7-meson have been challenged, and even in this case the situation remains rather 
controversial. The early experiments CELLO and CLEO 31], which studied the small-Q 2 behavior of the form factor 
-FM77, found it to be consistent with the asymptotic form, (f> as {z) = 6\/2/mz(1 — z). Later the BABAR collaboration [32j | 
found a rapid growth of the form factor Q 2 |-F,r 77 (Q 2 ) | 2 in the large-Q 2 regime. This observation drew attention to this 



problem and gave birth to speculations that the pion DA could be far from the asymptotic shape [51( (see also a recent 
review by Brodsky et. al. in [34], H^]). However, the most recent data from BELLE [33| did not confirm the rapid growth 
found by BABAR. As was found in [52|, [53[, the Gegenbauer expansion coefficients of the pion DA (t>2;ir{z) are small 
and at most give a 10% correction for the minus-first moment, based on the fits of BELLE, CLEO and CELLO data. 
For kaons there is no direct measurements of the DAs, it is expected however that the deviations from the pion DA are 
parametrically suppressed by the quark mass m s /GeV. Numerically this corresponds to a 10-20% deviation. 
For this reason in what follows we assume all the Goldstone DAs to have the asymptotic form, 

02;{ 7 rJC, ) }(z) ~ <f>as{z) = 6\/2/ M 2;(l - z) . 

For the decay couplings we use standard values ~ 93 MeV, Jk ~ 113 MeV, « fx. 

More than a dozen of different parametrizations for GPDs have been proposed so far [ll fl6[ ^1 59l| . While we neither 
endorse nor refute any of them, for the sake of concreteness we select the parametrization [54ll60l. l6l|. which succeeded to 
describe HERA [62| and JLAB H3, HH data on electro- and photoproduction of different mesons, so it might provide 
a reasonable description of z^DVMP. The parametrization is based on the Radyushkin's double distribution ansatz. It 
assumes additivity of the valence and sea parts of the GPDs, 



which are defined as 



n val 
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H q = 
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\a\ + \P\<l 



1 



.;ln|/3|)t. 



and q va i and q sea are the ordinary valence and sea components of PDFs. The coefficients bi, oti, as well as the parametriza- 
tion of the input PD Fs q ix), Aq(x) and pseudo-PDFs e(x), e(x) (which correspond to the forward limit of the GPDs E, E) 
are discussed in [54l. |6& [6lj . The unpolarized PDFs q(x) in the limited range Q 2 < 40 GeV 2 roughly coincide with the 
CTEQ PDFs. Notice that in this model the sea is flavor symmetric for asymptotically large Q 2 , 



(11) 
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FIG. 2: (color online) Pion production on nucleons. (a) Processes without strangeness production, (b) Processes with nucleon to 
hyperon transition (A5* = 1). Kinematics t — t m i n (Aj_ = 0) is assumed. 



where 

>^ = 1+ TToMmry Ql = iGeV2 - 

The equality of the sea components of the light quarks in (ITT|) should be considered only as a rough approximation, 
since in the forward limit the inequality d ^ u was firmly established by the E866/NuSea experiment [63| . For this reason, 
predictions done with this parametrization of GPDs for the p <^ n transitions in the region xsj € (0.1. ..0.3) might slightly 
underestimate the data. 



IV. NUMERICAL RESULTS AND DISCUSSION 



In this section we perform numerical calculations of the cross-sections of the processes listed in Table U relying on the 
GPDs presented in the previous section. The results for neutrino-production of pions on nucleons are depicted in Figure 

In the left pane of the Figure [2) we grouped the pion production processes without excitation of strangeness. The 
results for the diagonal channels, p — > p and n — > n extend to large xbj our previous calculations [64| for diffractive 
neutrino-production of pions performed in the dipole approach, which assumes dominance of the sea. Differently from 
sm&l\-x Bj diffraction, in the valence quark region we found that the production rate of 7r + on neutrons is about twice 
larger than on protons. This results from the fact that the handbag diagram in the proton probes the GPD H,i, whereas 
larger H u contributes via crossed handbag; in the case of neutron they get swapped. At large XBj <^ 0.6 the corresponding 
cross-section is suppressed due to increase of \t m in (x, Q 2 ) \ - 

The off-diagonal processes with p ^ n transitions are suppressed at small x Bj because they probe the GPD difference 
H u — Hd- In the small- a; regime (xsj ^0.1) the density of light sea quarks become equal, d ~ u, and cancel. The 
valence quark PDFs and the invariant amplitude Tm behave like ~ 1/^/XbJ, so that the cross-section vanishes as ~ xsj- 
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This result agrees with the Regge phenomenology, which predicts this cross-section to fall as 6 ~ x B , 2ap ^ XBj [49]. The 
cross-sections of the neutral ir° production on the proton and neutron (processes vp — > vir°p and vn — > vir°n) coincide 
under the assumption of iJ-dominance, however in the general case they differ, with effects ~ H u — Hd- Numerically 
in the considered parametrization of GPDs these effects are of order 1 %, so the difference between the two curves is 
invisible in the plot. A similar result holds for the processes vp — > vir + n and vn — > vn~p: the corresponding cross-sections 
exactly coincide under the assumption of ii-dominance but differ in the general case. As in neutral pion production, the 
difference is controlled by a small ~ H u — Hd, however, due to suppression of the GPD H in the small- x Bj region the 
effects proportional ~ H u — Hd are relatively large, and the difference between the two cross-sections becomes visible in 
the plot. 

In the right pane we show the cross-sections of pion production with nucleon to hyperon transition. These cross-sections 
are Cabibbo-suppressed and hardly can be detected in the Minerva experiment. In contrast to p n processes, at small 
XBj the sea contributes to the difference H u — H s , H d — H s . First of all, the sea flavor asymmetry appears due to the 
presence of proton nonperturbative Fock components, like p — > KA. This asymmetry vanishes in the invariant amplitude 
at small xsj as x B " K * , where the intercept of the K* Reggeon trajectory is ax* ~ 0.25. Correspondingly, this contribution 
to the cross section is suppressed as ~ x B j. 

At very small XBj a more important source of flavor asymmetry is the Pomeron itself. This is easy to understand within 
the dipole approach, in which the sea-quark PDF probed for instance by a virtual photon, corresponds to the transition 
7* — > qq and interaction of the qq dipole with the target proton (in the target rest frame). This cross section is different for 
ss and uu(dd) dipoles due to so called aligned-jet configurations, in which the dominant fraction of the dipole momentum 
is carried by the quarks. Such a non-perturbative contribution to the PDF is known to persist even at high Q 2 [65|. Such 

a flavor asymmetry of the sea rises in the cross section as x B j , controlled by the value of the Pomeron intercept, 
e (Q 2 ) = a P ~ 1) which was found in the experiments at HERA [66[ to grow steeply with Q 2 . 

At medium-small values of XBj in the transition region between the two above regimes the interference of these two 

mechanisms, which behaves as x\^ K * e ^ \ is also important. This region of xsj ^ 0.1 is the domain of our main 
interest, while small XBj *C 1 is beyond the scope of this paper. At XBj ^ 0.2, the dominant contribution comes from the 
valence quarks. 

In Figures 13141 we show the results for kaon production. Figure [3] presents the results for Cabibbo-allowed processes 
(AS* = 0), which involve a nucleon into the hyperon transition. For moderate values of xsj ^ 0.2 the dominant contribution 
to these processes comes from light quarks, and in the region of smaller xsj from s-quarks. In the left and right panes of 
Figure [31 we show the cross sections for neutral and charge kaon production respectively. The order of magnitude of these 
cross sections is comparable with that for pions, so should be easily measured at Minerva. However, all these processes are 
suppressed at high energies ( small- x regime). As was discussed in Section [Til there are two parametrization-independent 
relations between the kaon-hyperon production cross-sections, da up ^ fll - k+t,+ = 4 da ^+ K o^.o and da i)n ^ ll + KL gS - = 
dcrvn^u- k+ s° ■ I n Figure|3]one can see the cross-section without hyperon formation. Such diagrams are Cabibbo suppressed 
(have AS* = 1) and have too small cross-sections, hardly detectable in current experiments. 

Figure [5] demonstrates the results for ^-production, both with AS = and AS = 1. The cross-sections for all such 
processes at moderate values of XBj are found to be alike and to have similar dependences on xsj- Numerically, they 
differ at most by an order of magnitude. At small XBj however, only the neutral current cross-sections vp(n) — > vr\p(n) 
survive, which depends on xbj with the chosen parametrization of GPDs as ~ x b 2 - 

The ^-integrated cross section of pion production, calculated for diagonal transitions and plotted in the left pane of 
Figure [5] demonstrates the features similar to the forward cross section depicted in Figure [2] 

The cross section is steeply falling toward large xsj due to increasing |t m , n |, whereas at small xsj it behaves similarly 



6 The pion intercept o^-(O) ~ 0, so the pion exchange contribution is suppressed as x\ 



9 



to 
CM 

J 10- 2 



10" 3 



> 10" 4 



y=0.5, Q 2 =4 GeV 2 



- 1 1 1 — r~ 



(a) 



0.1 



vp^n + K° LS i: 

vp^(i + K° L ' s A 

vn^(x + K°' L s £- 
viwvK° L S A 
vn^vK° L ' s i: 




0.2 0.3 0.4 0.6 0.8 



10 -1 

to 

55 

CM 

J 10- 2 



10-3 



> 1Q-4 



~l 1 1 1 1 1 1 1 1 - 

y=0.5, Q 2 =4 GeV 2 (b) 



vp^n"K + i: + 

vp->vK + A 

vp^vK + £° 

vn^n-K + E° 

vn->|i"K + A 

vn-wK+r 




0.1 



0.2 0.3 0.4 0.6 0.8 



FIG. 3: (color online) Kaon production on the nucleon with nucleon to hyperon transition (AS = 0). (a) Neutral kaon production, 
(b) Charged kaon production. Kinematics t = t m in (A_l = 0) is assumed. 
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FIG. 4: (color online) Kaon production on the nucleon without nucleon to hyperon transition (AS = 1). Kinematics t — tmin 
(Ax = 0) is assumed. 



to the unintegrated cross section. Although in the valence region the cross-sections may differ up to a factor of two, all 
diagonal channels for charged and neutral currents unify at the same production rate at small xsj , confirming the results 
of the dipole description [64| . 

The t-dependence of the differential cross-section is controlled by the underlying parametrization of GPDs we rely upon, 
and our results for the differential cross section of neutrino-production of pions arc plotted in the right pane of Figure [S] 
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FIG. 5: (color online) ^-production on the nucleon. In the upper right corner are the Cabibbo-allowed processes (AS = 0), in the 
lower left corner are Cabibbo-suppressed processes (AS = 1). Kinematics t = t m i n (Aj_ = 0) is assumed. 
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FIG. 6: (color online) (a) f-integrated two-fold cross-section d?a/dvdQ 2 . (b) f-dependence of the differential cross-sections for 
selected processes. 



It can be roughly approximated by the exponential t-dependence da ~ exp(B e fft), with the slope 

1 duyyN—^MN 



B eff = 



aw 



dt 



which decreases with xbj from about 6 down to 2 GeV 2 , as is shown in the left pane of Figure [7] To a good extent the 
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x b\ Q 2 , GeV 2 

FIG. 7: (color online) Effective slope -B e // which controls the i-dependence of the differential cross-section, as a function of XBj 
(left) and virtuality Q 2 (right). 

calculated xsj-dependence of the slope is described by B(x) = Bq — /3ln(x), where the coefficient j3 may vary between 0.3 
(« 2a' sea (0)) and 1.8 (= 2a' val (0)), depending on the process and value of Q 2 . 

The Q 2 dependence of B e ff depicted in the right pane of Figure [7] is rather mild, which is due to the weak dependence 
of the shape of PDFs on Q 2 under DGLAP evolution. 

Finally, in order to demonstrate explicitly the effect of skewness we compare in Figure [5] the cross-sections of several 
processes calculated in the model [5J, l60l Ifillj , and with the simple zero-skewness parametrization 

H f (x,Z,t)f*q f (x)F N (t), (12) 

where qf(x) is the parton distribution, and F^(t) is the nucleon form factor. We see that the results of the two parametriza- 
tions differ up to a factor of two. 

V. SUMMARY 

We evaluated the cross-sections of deeply virtual meson production for pions, kaons and eta-mesons in neutrino-nucleon 
interactions. The production rate of the Cabibbo-allowed processes in the Bjorken regime is found to be sufficiently 
large to be detected at the current level of statistics of neutrino experiments, in particular in the Minerva experiment at 
Fermilab, with accuracy, which allows to disentangle between different models of GPDs. For this purpose we provided 
detailed information on the distributions of the production rate versus different variables, xsj t, Q 2 . For further practical 
applications we provide a computational code. We also evaluated the cross sections of Cabibbo-suppressed channels 
(AS = 1), but found them too weak to be detected by any of forthcoming experiments. 
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FIG. 8: (color online) vDVMP cross-section for certain processes. The abbreviations KG and ZS stand for the Kroll-Goloskokov 



After this manuscript was submitted, we learned that neutrino-production of pions was also studied in 67] within 
the approach proposed in [68j . The central assumption of [67j is that the pion coupling has a q^q structure, i.e. the 
dominant contribution comes from the subleading twist pion DA <f> p . This contribution represents a 0(1/Q 2 ) correction 
to the longitudinally polarized cross-sections discussed in this paper, however it may be numerically important for the 
transversely polarized current. Then, the amplitude gets contributions from transversity GPDs Ht, Et, Ht, Et- We 
did not include such corrections, because apart from the uncertainty in the chiral odd GPDs, this requires modelling of 
the poorly known twist-three pion DA 4> p . No model-independent estimate for this quantity is available so far, while 
mo del- dependent results vary considerably (69l . l70T | . 
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